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1  Introduction 


In  the  life  of  a  eombat  soldier,  traumatic  cataract  in  ocular  lenses  may  result  from  blast  load¬ 
ing,  whereby  (i)  the  lens  capsule  (Fig.l)  is  perforated  by  intraocular  foreign  bodies  (lOFBs  [Wal¬ 
ter,  1962,  Mader  et  ah,  1993,  Parver  et  ah,  1993,  Wong  et  ah,  1997,  Mader  et  ah,  2006,  We- 
ichel  and  Colyer,  2008])  which  in  turn  damage  the  lens  fiber  cells,  (ii)  the  lens  is  loaded  fluid 
dynamically  by  the  surrounding  aqueous  and  vitreous  humors  [Banitt  et  ah,  2009]  (see  Fig.l), 
and/or  (iii)  the  lens  internal  substance  (crystallins  lens  fiber  cells)  is  stressed  by  the  passing 
shock  wave.  Traumatic  cataract  can  result  in  a  partially  or  fully  clouded  lens,  complete  dislo¬ 
cation  of  the  lens  (floating  between  aqueous  and  vitreous  humors,  see  Fig.l),  or  zonule  rupture 
such  that  partial  or  full  vision  loss  may  occur.  The  mechanisms  of  traumatic  cataract  forma¬ 
tion  that  may  require  cataract  surgery  (implantation  of  an  intraocular  lens  (lOL))  are  not  well 
understood  in  comparison  to  the  mature  and  ever-improving  surgical  technology  and  procedures. 
The  hypothesis  of  the  research  is  that  an  ultrastructurally- 
based  computational  finite  element  model  of  the  ocular 
lens  subjected  to  blast  loading  can  assist  in  better  under¬ 
standing  how  traumatic  cataract  is  formed  in  the  combat 
soldier,  and  in  turn  improve  our  understanding  of  traumatic 
cataract  in  civilians  whose  eyes  are  subjected  to  impact 
loading.  The  scope  of  the  research  is  to  develop  a  multi¬ 
scale,  ultrastructurally-based,  computational  model  of  the 
ocular  lens  subjected  to  blast  loading,  in  conjunction  with 
imaging  methods  to  identify  lens  capsule  and  internal  sub¬ 
stance  structure  and  mechanical  experiments  for  calibrat¬ 
ing  material  model  parameters. 


Figure  1.  Eye  cross-section. 
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2  Body 

The  research  tasks  focussed  on  over  the  past  year  have  been  the  imaging  ones  (Task  6),  and  some 
modeling  (Task  1-3).  For  Task  6,  undergraduate  students  Sai  and  Sri  Radhakrishnan  have  taken 
the  lead  on  the  cryo-electron  microscopy  and  tomography  (cryo-EM  and  ET)  (Sai)  imaging  of  type 
IV  collagen  structure,  and  confocal  laser  scanning  microscopy  (CESM)  (Sri)  for  the  internal  lens 
fiber  cells.  In  addition,  on  the  modeling  side  a  PhD  student.  Boning  Zhang,  has  been  recruited  to 
work  on  Tasks  1-3.  New  results  for  a  combined  Eulerian-Eagrangian  approach  to  large  deforma- 
tional  loading  of  the  lens,  and  solid-fluid  interaction,  are  presented  in  collaboration  with  Assistant 
Professor  Eranck  Vernerey  and  his  PhD  student,  Eouis  Eoucard. 

A  No  Cost  Extension  (NCE)  will  be  filed  in  May  2014,  to  extend  the  project  through  September 
2016. 

Next,  research  on  the  Tasks  over  the  past  year  are  summarized. 
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6a.  Before  and  after  whole  porcine  lens  unconfined  compression,  image  lens  fiber  cell  geom¬ 
etry  using  CLSM  (low  strain  rate  subtask  4a:  months  1-4;  higher  strain  rate  sub  task  4b: 
months  6-8). 

In  collaboration  with  Dr.  Christopher  English  of  the  Moleeular,  Cellular,  and  Developmental 
Biology  (MCDB)  department,  we  are  still  refining  the  procedure  of  CLSM  for  identifying  the 
ultrastructure,  at  various  cross-sections  through  the  whole  lens,  of  the  internal  lens  fiber  cells.  We 
are  foeussing  our  efforts  on  as-reeeived  poreine  lenses  and  one  human  lens. 

We  have  eondueted  preliminary  analysis  on  confocal  images  of  human  interior  lens  substance, 
attempting  to  identify  the  ultrastructure  of  the  lens  fiber  eells.  The  first  image  of  the  sliee,  and 
the  location  of  the  sliee  in  the  eontext  of  a  drawing  of  the  lens  fiber  eell  ultrastructure  (drawing 
from  pg99  of  [Kessel  and  Kardon,  1979])  is  shown  in  Fig. 2.  The  donor  was  a  66-year-old  male, 
with  lens  obtained  in  Balaneed  Salt  Solution  (BSS)  within  48hrs  post  mortem.  When  the  lens  was 
prepared  for  eonfoeal  laser  seanning  mieroseropy  (refer  to  Seetion  7.1  of  Christopher  Bays  MS 
thesis,  ineluded  as  an  attaehment  to  Year  2  Final  Report),  the  nueleus  and  eortex  separated,  and 
thus  in  the  image  only  the  eortex  seetion  is  shown. 


Phalloidin  Stain-  actin 


(fiber  cells) 


Figure  2.  CLSM  slice  of  human  lens  fiber  cell  structure  aligned  with  approximate  location  within  full  lens 
(drawing  from  pg99  of  [Kessel  and  Kardon,  1979]). 
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Figure  3  shows  zoomed-in  images  of  the  lens  fiber  eell  ultrastrueture.  When  eompared  to 
images  for  the  2+yr  old  pig  lens  (see  Chapter  7  of  Christopher  Bays  MS  thesis,  as  attaehment  to 
the  Year  2  Final  Report),  the  66-year-old  male  fiber  eell  ultrastrueture  is  less  identifiable.  One 
hypothesis  for  this  less-identifiable  fiber  cell  ultrastrueture  in  the  older  human  lens  is  the  wearing 
or  degradation  of  the  lens  fiber  cells  over  longer  time.  This  hypothesis  remains  to  be  tested  as  more 
data  are  obtained  for  younger  human  lenses,  and  from  other  older  human  lenses. 


Figure  1:  lOX  confocal 
image  of  human  lens 
(image  superimposed  in 
lens  model) 

Figure  2:  20X  (region  of 
human  lens  outlined  by 
red  square) 

Figure  3:  20X  (region  of 
human  lens  outlined  by 
purple  square) 

Figure  4:  20X  (region  of 
human  lens  above  that 
shown  in  lOX  image) 


Figure  3.  Zoomed-in  images  of  CLSM  slice  of  human  lens  fiber  cell  structure. 
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6b.  On  as-received  porcine  lens  capsules,  image  type  IV  collagen  ultrastructure  in  lens  cap¬ 
sule  using  cyro-electron  tomography  (months  1-4). 

Identification  of  type  IV  collagen  meshwork  ultrastructure  in  the  lens  capsule  using  cyro- 
electron  microscopy/tomography  is  described  in  Chapter  6  of  C.  Bay’s  thesis,  included  as  Ap¬ 
pendix  to  last  year’s  annual  report.  The  current  focus  of  the  research  is  to  determine  if  the  images 
we  are  obtaining,  such  as  in  Fig. 4,  are  indeed  structure,  or  if  they  are  an  artifact  of  the  freezing  pro¬ 
cess.  We  are  collaborating  with  Dr.  Thomas  Giddings  and  Dr.  Andreas  Hoenger  of  the  Molecular, 
Cellular,  and  Development  Biology  (MCDB)  Department  at  the  University  of  Colorado,  Boulder, 
about  this  topic.  More  experiments  are  being  conducted  to  determine  whether  what  we  are  seeing 
is  structure  or  artifact,  such  as  increasing  the  percentage  of  cryo-protectant  solution. 


capsule-Epon-lSOmM  mannitol  Epon-mesh-3m.tif 
274/SE  Epon  block  F 
capsule-  mesh?  8  nm  fibers 
Print  Mag:  223000x  @  7.0  in 

1:  31:  49  PM  9/10/2013 


100  nm 
HV®100 . OkV 
Direct  Mag:  ISOOOOx 
AMT  Camera  System 


Figure  4.  Cryo  EM  image  of  what  is  possibly  type  IV  collagen  structure. 
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2c:  Formulation  and  finite  element  implementation  of  multiscale  perforating  finite  strain 
biphasic  mixture  (solid  and  fluid)  solid-shell  continuum  model  of  lens  capsule  in  Tahoe, 
to  model  cutting/  perforation  of  lens  capsule  based  on  implementation  in  subtasks  Ic  and 
2b  (months  28-36).  3b:  Using  result  of  subtask  3a,  formulate  and  implement  multiscale 
hierarchical,  anisotropic,  lens  fiber  cell  equivalent  soft  viscoelastic  constitutive  model  of  the 
internal  lens  substance  (months  40-  48). 

Large  Deformation  Lagrangian-Eulerian  Computational  Modeling  of  Lens  Puncture  and  Also 
Solid-Fluid  interaction  between  Lens  and  Surrounding  Fluids: 

For  computational  modeling,  we  are  attempting  a  large  deformation  hybrid  Lagrangian- 
Eulerian  simulation  of  lens  puneture  tests  in  order  to  eventually  model  penetration  by  Intra-Ocular 
Foreign  Bodies  (lOFBs)  (Seeondary  Blast),  but  also  shoek  propagation  and  solid-fluid  interaetion 
between  the  lens  and  vitreous  and  aqueous  humors  (Primary  Blast).  This  is  in  eollaboration  with 
Assistant  Professor  Franek  Vemerey,  and  his  PhD  student  Mr.  Louis  Foueard,  both  at  CU  Boulder. 
The  images  below  illustrate  a  preliminary  simulation  of  spherieal  tip  indentation  into  a  2-i-yr  old 
pig  lens. 

The  method  takes  an  Eulerian  approaeh  to  deseribe  the  large  deformations  of  an  elastic  mem¬ 
brane  and  its  interactions  with  the  surrounding  fluid.  The  membrane  is  modeled  as  a  two- 
dimensional  elastic  surface,  aeross  whieh  diseontinuities  of  continuum  fields  sueh  as  pressure  and 
fluid  veloeity  in  the  tangential  direction  can  naturally  be  enforeed  using  the  extended  finite  element 
method  (X-FEM).  The  tracking  and  evolution  of  the  membrane  is  handled  with  the  Grid  Based  Par- 
tiele  method  that  is  well  suited  to  evaluate  the  membrane  higher  order  geometrieal  information. 

The  results  of  the  simulation  of  indentation  of  a  pig  lens  by  a  spherieal  indenter  are  eompared 
with  experimental  data.  The  experimental  data  are  obtained  for  a  2-i-  year-old  poreine  lens  under¬ 
going  indentation  by  a  2mm  (in  diameter)  indenter  with  a  spherieal  tip  (Pig. 5).  The  indenter  is 
lowered  at  a  speed  of  0.5mm/s  until  failure/puneture  of  the  eapsule.  The  simulation  was  earried 
out  using  a  newly  developed  Eulerian  approaeh  able  to  deseribe  the  large  deformations  of  an  elas- 
tie  membrane  and  its  interaetions  with  the  surrounding  and  enelosed  fluid  (Pig. 5).  The  eapsule 
membrane  is  assumed  to  be  ineompressible  and  isotropie,  for  the  time  being.  The  Mooney-Rivlin 
strain-energy  funetion  was  chosen  to  model  membrane  elastie  response,  in  Eq.(l) 

^  ^  2(1  +  a)  ~  -1- -h  Aj  ^^2  ^  —  3)  H-  a(Xf  -l- -l- Af  A|  —  3)]  (1) 

where  /i  is  the  infinitesimal  shear  modulus  of  the  elastic  material,  Ct  is  a  parameter  that  controls  the 
deviation  from  a  linear  Hookean  response,  and  Ai  and  A2  are  the  prineipal  stretehes  in  direction  1 
(radial)  and  2  (axial) .  In  the  simulation,  the  eapsule  is  eonsidered  filled  with  a  Newtonian  viseous 
fluid  of  viscosity  /if. 

In  this  simulation,  the  indenter  (blaek  cylinder)  in  Pig. 5  progressively  approaches  and  deforms 
the  elastie  eapsule  filled  with  a  viseous  fluid  and  resting  on  a  plane  surface.  The  initial  geometry  of 
the  eye  lens  was  extraeted  from  the  experimental  video  of  indentation  at  time  t  =  0,  and  the  system 
was  meshed  aeeordingly  using  a  23x30  element  grid.  The  (frietion-less)  interaetion  between  the 
eapsule  and  the  indenter  is  handled  via  the  penalty  method.  Pigure  6  shows  the  dilatation  (jaeobian 
J)  of  the  elastie  membrane.  Pigure  6  shows  the  deformation  gradient  eomponents  in  the  axial  (F22) 
and  radial  (Pn)  direetions.  One  ean  observe  the  very  large  deformation  (approaching  800%)  and 
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Figure  5.  Indentation  and  large  deformations  of  the  porcine  (2+years  old)  lens:  (left)  experiment,  (right) 
simulation.  The  indenter  has  a  2mm  diameter  and  a  spherical  tip  and  is  lowered  at  0.5mm/s  until  puncture 
of  the  capsule. 


dilatation  of  the  membrane  at  the  tip  of  the  indenter,  whieh  would  prove  very  diffieult  to  simulate 
using  a  Lagrangian  approaeh.  This  stress  eoneentration  at  the  tip  of  the  indenter  is  also  eonsistent 
with  the  faet  that  the  rupture  of  the  eapsule  membrane  was  observed  to  take  plaee  at  the  same 
loeation  experimentally.  Finally,  Fig. 6  shows  the  normalized  foree  versus  displaeement  eurve  for 
both  the  simulation  and  the  experimental  results.  Up  to  the  point  of  failure  of  the  eapsule  (peak), 
the  simulation  results  show  good  eonsisteney  with  the  experimental  data. 
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Force  VS  displacement 


Figure  6.  Indentation  and  large  deformations  of  lens,  (top  left)  dilatation  of  the  surface  (Jacobian),  while 
(top  right)  and  (bottom  left)  show  deformation  gradient  in  axial  {F22)  and  radial  (Fn)  directions,  (bottom 
right)  normalized  force  versus  displacement  curve  for  both  the  simulation  results  and  the  experimental  data. 
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3  Key  Research  Accomplishments 


•  CLSM  images  of  the  internal  lens  fiber  eell  strueture  are  showing  promise,  and  we  are  eonfi- 
dent  we  ean  obtain  good  struetural  identifieation  of  lens  fiber  eell  strueture  pre  and  post-test. 
This  will  be  used  to  identify  the  lens  fiber  eell  strueture  for  the  multiseale  modeling  of  the 
internal  lens  substanee. 

•  Eleetron  Mieroseopy  (EM)  and  Tomography  (ET)  images  are  not  yet  believed  to  be  strueture, 
but  eould  possibly  be  an  artifaet  of  the  freezing  proeess.  This  is  beginning  to  raise  questions 
as  to  the  validity  of  type  IV  eollagen  struetural  images  provided  in  the  literature  to  date.  The 
MCDB  department  and  Bio3D  lab  bio3d.colorado.edu  are  world-renown,  NIH-funded 
faeilities  for  identifying  strueture  of  biologieal  tissues.  Thus,  they  have  exeellent  faeilities, 
and  are  eareful  in  their  tissue  preparation  and  interpretation  of  results. 

•  Preliminary  large  deformation  Eagrangian-Eulerian  axisymmetrie  eomputational  modeling 
of  lens  indentation,  and  also  solid-fluid  interaetion  between  the  lens  and  surrounding  fluids, 
shows  promise  for  providing  a  robust  large  deformation  eomputational  meehanies  frame¬ 
work  for  simulating  Primary  and  Seeondary  blast  on  the  lens,  as  well  as  the  whole  eye 
globe,  potentially. 

4  Reportable  Outcomes 

1 .  August  2013,  Experiments  and  Einite  Element  Analysis  of  Stress  Relaxation  upon  Compres¬ 
sion  of  Whole  Poreine  Eenses,  EMI  2013,  Northwestern  University. 

2.  On  June  24-25,  Dr.  Matthew  Reilly  from  the  University  of  Texas,  San  Antonio  (UTS A), 
visited  PI  Regueiro  to  diseuss  eaeh  other’s  researeh  and  how  we  ean  work  together. 

3.  The  eollaboration  beginning  with  Dr.  Vernerey  of  CU  Boulder  on  large  deformation 
Eagrangian-Eulerian  eomputational  modeling  of  lens  indentation,  and  also  solid-fluid  inter- 
aetion  between  the  lens  and  surrounding  fluids,  shows  promise  for  providing  a  robust  large 
deformation  eomputational  meehanies  framework  for  simulating  Primary  and  Seeondary 
blast  on  the  lens,  as  well  as  the  whole  eye  globe,  potentially. 

5  Conclusion 

The  researeh  progress  over  the  last  year  has  foeussed  on  imaging  of  the  lens  fiber  eell  strueture 
using  eonfoeal  laser  seanning  mieroseopy  (CESM),  and  type  IV  eollagen  strueture  using  eryo- 
eleetron  mieroseopy  and  tomography  (eryo-EM  and  ET).  The  CESM  promises  to  generate  good 
struetural  identifieation  of  lens  fiber  eells,  while  the  eryo-EM  and  ET  have  not  yet  positively  iden¬ 
tified  type  IV  eollagen  strueture  of  the  lens  eapsule. 

The  eollaboration  with  Dr.  Vernerey  of  CU  Boulder  on  large  deformation  Eagrangian- 
Eulerian  eomputational  modeling  of  the  lens  will  be  extended  to  three-dimensions  (3D)  and  ultra- 
strueturally-based  modeling  of  eapsule  and  fiber  eells.  Sueh  modeling  eould  be  extended  to  the 
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whole  eye  globe,  given  the  ease  of  using  Eulerian-based  methods  in  large  deformation  eomputa- 
tional  meehanics  and  solid-fluid  interaetion.  Traeking  interfaees  is  the  ehallenge,  and  sueh  details 
are  being  investigated  as  part  of  the  researeh. 
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